Fermionic quantum dimer and fully-packed loop models on the square lattice 



O 



Frank Pollmann,' Joseph J. Betouras,^ Kirill Shtengel,^ and Peter Fulde' '* 

' Max-Planck-Institute for the Physics of Complex Systems, D-01187 Dresden, Germany 
^Department of Physics, Loughborough University, Loughborough, Leicestershire LEll 3TU, UK 
'^Department of Physics and Astronomy, University of California, Riverside, CA 9252L USA 
^The Asia Pacific Center for Theoretical Physics, Pohang Korea 
(Dated: January 4, 201 1) 

We consider fermionic fully-packed loop and quantum dimer models which serve as effective low-energy 
models for strongly correlated fermions on a checkerboard lattice at half and quarter filling, respectively. We 
identify a large number of fluctuationless states specific to each case, due to the fermionic statistics. We discuss 
the symmetries and conserved quantities of the system and show that for a class of fluctuating states in the half- 
filling case, the fermionic sign problem can be gauged away. This claim is supported by numerical evaluation of 
the low-lying states and can be understood by means of an algebraic construction. The elimination of the sign 
problem then allows us to analyze excitations at the Rokhsar-Kivelson point of the models using the relation to 
the height model and its excitations, within the single-mode approximation. We then discuss a mapping to a U( 1 ) 
lattice gauge theory which relates the considered low-energy model to the compact quantum electrodynamics in 
2-1-1 dimensions. Furthermore, we point out consequences and open questions in the light of these results. 
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I. INTRODUCTION 

Strongly correlated electrons on lattices with geometric 
frustration can lead to new physical properties. These come 
essentially from the quantum fluctuations acting on an exten- 
sively degenerate collection of classical ground states. From 
that point of view, strongly interacting spinless fermions (or 
alternatively, fully spin polarized electrons) on a frustrated lat- 
tice present an interesting example of this kind of physics. A 
spectacular consequence for example is that at certain fillings, 
fractionally charged excitations have been predicted for the 
case of large nearest-neighbor repulsion on certain frustrated 
lattices J. The subject is not purely academic since there is ex- 
perimental evidence that electrons in systems which have the 
considered lattice structure can be strongly correlatedi^^ Fur- 
thermore, certain optical lattices with frustrated interactions 
can be produced by interfering laser beams and loaded with 
fermionic atoms. By varying the lattice constant, the strength 
of the correlations of the particles can be modified^i 

In the present work we focus on the problem of strongly 
correlated electrons on a checkerboard lattice at two different 
filling fractions. This lattice can be seen as a two-dimensional 
projection of the pyrochlore lattice (3D tetrahedra are repre- 
sented by 2D squares whose diagonals and sides are on an 
equal footing as far as couplings are concerned). We will 
consider the case of strong nearest neighbor repulsion that 
favors the least possible number of fermions in each such 
square. At half-filling, this number is two, and an effective 
low-energy model maps onto a quantum fully-packed loop 
(FPL) model on a square lattice. Some parts of this work 
have been briefly discussed before;^ In the case of quarter- 
filling, a lowest potential energy arrangement corresponds to 
one fermion per square which can be mapped onto a hard- 
core dimer model. Here we supplement the study with fur- 
ther results and different aspects of the problem regarding the 
symmetries and conserved quantities, the sign problem and 
the corresponding gauge theories. We will discuss a num- 



ber of consequences, especially related to the fermionic na- 
ture of the loops/dimers in this particular problem. Quantum 
dimer models (QDMs) have have been studied quite exten- 
sively in the past two decadesi^ Originally introduced in the 
context of quantum magnetisn^, they became a focus of re- 
newed interest following the discovery of a gapped quantum 
liquid phase on a triangular lattice.^ It has been further estab- 
lished that a gapped topological phase with deconfined exci- 
tations generically exists in 2D for QDMs on non-bipartite 
latticesi^— Such a phase has been found to have an effective 
description in terms of a Z2 gauge theory»i^iii By contrast, 
on bipartite lattices (e.g. a square lattice) there is a single 
deconfined quantum critical point separating phases of bro- 
ken translation/rotation symmetry i^ii^^ Such a point is char- 
acterized by gapless excitation spectrum which in turn leads 
to logarithmic confinement at any finite temperature. In both 
cases, QDMs can be mapped onto loop models, with one sig- 
nificant difference: on a bipartite lattice the resulting loops 
can be naturally oriented while there is no natural notion of 
direction on a non-bipartite lattice. As a consequence, an ef- 
fective gauge theory for a bipartite case is now a U(l) the- 
ory (the oriented loops can be though of as field lines in the 
absence of sources)fi^iii In addition, oriented loops naturally 
lend themselves to a height representation whereby they serve 
as contour lines. Several related models, such as the quantum 
six/eight-vertex models in 2D have been shown to conform to 
the same dichotomy: a model with orientable loops is either 
ordered or critical and is described by a U(l) gauge theory 
while non-orientable quantum loops can result in a decon- 
fined topological phaseji^— In all these models, the Hilbert 
space typically separates into disconnected subspaces corre- 
sponding to various conserved quantum numbers. When these 
quantum numbers are of non-local topological nature, such as 
winding numbers, the corresponding subspaces are refeiTed to 
as topological sectors. While a Z2 topological phase is char- 
acterized by four such sectors and a corresponding four-fold 
ground state degeneracy on a torus, the U(l) models have a 
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much larger number of conserved quantities. 

All aforementioned models share one important feature: 
all off-diagonal matrix elements connecting various quantum 
dimer/quantum loop states are non-positive (or can be made 
non-positive by a trivial gauge transformation such as associ- 
ating a phase to dimers at certain locations). Within each dis- 
connected sector the Hamiltonian dynamics is ergodic in the 
following sense: every two quantum states form a non-zero 
matrix element with some power of the Hamiltonian. (Here 
we are not concerned with the issue of how such a power 
may depend on the system size - a question essential for un- 
derstanding the spectral gap or a lack thereof.) The Perron- 
Frobenius theorem can then be applied to exp(— ri/ ) for some 
positive T within each disconnected subspace; it is follows 
that the ground state for each sector is unique and nodeless.i^i^^ 

Much less is known about models with non-Frobenius dy- 
namics which results in some positive matrix elements. A 
striking difference between Frobenius vs. non-Frobenius 
QDMs on the kagome lattice was found in Refs. 1 12.23ll: 
while the conventional QDM exhibits a gapped Z2 topologi- 
cal phase, its counterpart with different signs has an extensive 
ground state degeneracy. A different class of non-Frobenius 
lattice models has also been studied recently: fermionic mod- 
els with supersymmetric Hamiltonians; in a number of cases 
such systems have been shown to have an extensive GS 
degeneracy.-^—. From these examples, it appears that the 
non-Frobenius dynamics are often the source of additional 
quantum frustration resulting in the extensive ground state en- 
tropy. 

The goal of the present manuscript is to present a number 
of results (and some open questions) related to another non- 
Frobenius model that was originally introduced in Refs. 
I28I1 . The paper is organized as follows: We begin by introduc- 
ing the considered model and derive an effective low-energy 
Hamiltonian in the following section. In Section|IIIl we begin 
by analyzing the conservation laws specific to both quarter 
and half-filled cases, and the constraints they impose on the 
quantum dynamics. We discuss the nature of the sign prob- 
lem and present two different approaches to alleviating it. In 
Section |IV] we extend the model by adding an extra term in 
the effective Hamiltonian that punishes flippable configura- 
tions. In such an extended model, we analyze in detail the 
Rokshar-Kivelson (RK) point and discuss the phase diagram. 
In Section|Vl we discuss the gauge theory for both the quarter- 
filled and the half-filled cases. We conclude and summarize 
our results in SectionlVll 

II. THE MODEL 

Our starting Hamiltonian for spinless fermions is of the fol- 
lowing form: 

(* 3) (i j) 

where the summation is performed over neighboring sites of 
a checkerboard lattice. We are interested in the strongly inter- 
acting limit V t > Q. Note that each site has six nearest 



neighbors as in the pyrochlore lattice. The configurations fa- 
vored by the interaction term depend on the doping, i.e., on 
the ratio of the number of particles to the number of the avail- 
able sites. The fact that the particles are spinless fermions 
automatically prohibits double occupancy of sites. Further- 
more, the interaction term favors uniform minimization of the 
number of such fermions on every crisscrossed square, a "pla- 
nar tetrahedron." We will be mostly interested in a case of 
half-filling whereby the potential energy term is minimized 
whenever there are two fermions on each "planar tetrahedron" 
as shown in Figure [lja,b) (the tetrahedron rulei^.). The other 
case considered in parallel with the half-filling in this work is 
that of quarter-filling whereby the interaction energy is min- 
imized by having one fermion per tetrahedron as shown in 
Figure [nc,d).= In either case, configurations satisfying these 
rules can be represented by FPLs / dimers on a square lattice 
connecting the centers of the planar tetrahedra. The particles 
are sitting in the middle of the bonds of the square lattice. At 
half-filling, the tetrahedron rule translates into a FPL cover- 
ing. By the same token, at quarter-filling, the minimum inter- 
action energy manifold is spanned by the dimer coverings of 
the square lattice such that each site belongs to one and only 
one dimer. In either case, all other configurations are sepa- 
rated by a potential energy gap of at least V . Different FPL 
/ dimer configurations are orthogonal to each other since they 
correspond to different locations of fermions (here we assume 
that these locations simply label the orthonormal set of Wan- 
nier functions). 



(a)Two-particle exchange at 
half-filling 




(c)Two-particle exchange at 
quarter-filling 



(b)Three-particle exchange at 
half-filling 




(d)Thi'ee-particle exchange at 
quarter-filling 



FIG. 1 : (a,b) - A configuration of fermions satisfying the "two per 
tetrahedron" rule at half-filling and its representation in terms of 
FPLs. (c,d) - Configurations of fermions satisfying the "one per 
tetrahedron" rule at quarter-filling and their representations in terms 
of quantum dimers. In both cases, "flippable plaquettes" correspond- 
ing to two- and three-particle ring exchanges are shaded (see text for 
details). 
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Both at half- and at quarter-filhng, the kinetic term in 
the Hamiltonian ([T]) moves particles between the neighboring 
sites thus creating configurations in which the respective tetra- 
hedron rules are violated. However, we are interested in the 
physics at the energy scale much smaller than V and therefore 
we should look at the correlated multi-particle hops (ring ex- 
changes) instead. The lowest-order process that does not take 
the system out of the low-energy manifold is a two-particle 
ring exchange depicted in Figure[na,c). Such a process repre- 
sents the dominant quantum dynamics for bosonic systems as 
well as for resonating valence bonds in spin systems.*^ How- 
ever, its amplitude vanishes identically for fermions: the rel- 
ative signs of clockwise and counterclockwise processes are 
opposite as a consequence of the fermionic exchange statis- 
tics, and the two contributions cancel each other Moreover, 
the same observation holds for all ring exchanges involving 
any even numbers of fermions; the contributions from clock- 
and counterclockwise moves of an even number of fermions 
similarly cancel. Therefore, the shortest allowed fermionic 
ring exchanges involve three particles rather than two (see 
Figure[ljb,d)).— Note that if the spin of fermions is included, 
processes of order t"^ /V do no longer vanish: two electrons 
with opposite spin can resonate around empty plaquettes of 
the checkerboard lattice and lower the energyj^— By project- 
ing the full Hamiltonian of Eq. ^ onto the Hilbert subspace 
restricted by the tetrahedron rules, we can perturbatively ob- 
tain an effective ring-exchange Hamiltonian that acts within 
such a subspaceiS To lowest non-vanishing order in t/V, such 
effective Hamiltonian becomes; 

E (I I +H.C.) 

{m,a} 

= 5E(l^)(^l + l^)(^l) (2) 
p 

at quarter-filling and 

= E (I :::n+H.c.) 

{m.B} 

-g E (|3:i)(i:C|+H.c.) 

{m,a} 

p 

at half-filling. Here the sums are performed over all polygons 
of perimeter six and g = 12t^/V^. As we have already men- 
tioned, the actual fermions are located at the centers of the 
dimers. The different relative signs of two terms at a half- 
filling signs result from the number of permuted fermions. It 
is either even or odd which in turn depends on whether the 
middle site of the hexagon is occupied or empty. By contract, 
all ring exchanges come with the same sign in the case of the 
quarter-filling, as they all correspond to an empty middle site. 

We should also point out that while the sign of the coupling 
constant g is determined by the sign of t, it can be effectively 
changed by the following local gauge transformation: mul- 
tiply a states by a factor of (—1) for each vertical dimer in 



columns numbered 4k and 4fc + 1 (where k G Z); do the 
same for horizontal dimers replacing columns by rows. A 
simple check shows that all ring exchange terms in Eqs. Q 
and ^ change their sign as a result. (It is worth noting that 
any diagonal, i.e. potential energy terms that one might con- 
sider adding to these Hamiltonians would remain intact under 
this gauge transformation.) Therefore, despite all appearances 
to the contrary, the quantum dynamics in Eq. Q is not non- 
Frobebious. Such a gauge transformation, however, does not 
fix the relative sign of the different ring exchange terms in 
Eqs. ^ - this issue will be addressed in the following sec- 
tion. 



III. QUANTUM DYNAMICS 

In this section, we discuss several aspects of the effective 
low-energy Hamiltonians which were derived in the previous 
section. 



A. Conserved quantities 

Our first observation is that both at quarter and half-filling, 
the low energy Hilbert space of our model can be represented 
by height configurations. This follows from the aforemen- 
tioned geometric representations of the ground states of the 
potential energy term in Eq. ([T]): in the case of quarter-filling 
they are the dimer coverings of the square lattice while for the 
case of half-filling they form FPL on the same lattice, which 
are isomorphic to six-vertex configurations Both, dimer 
coverings and six-vertex configurations can in turn be mapped 
onto configurations of an interface - a so-called height repre- 
sentation, also known as a solid-on-solid (SOS) model.—— 
Figure|2ja,b) show one of the ways of assigning heights to the 
plaquettes of the square lattice for both dimer coverings and 
FPL configurations. 

The quantum dynamics that preserves these geometric con- 
straints then results in quantum fluctuations of the local height 
variable. Hence our model is a variant of a quantum height 
model.i^^^ Let us begin by reiterating some results known 
from studies of bosonic quantum dimer models^Si^ and quan- 
tum six-vertex models i'^-'^'^^ A conservation law which is 
valid for the ring exchanges of any finite length is the conser- 
vation of the global slope of the height field. The global slope 
variable has two components, Kx{y) and Ky{x) in x and y di- 
rections, which form a set of conserved numbers. (In case of 
periodic boundary conditions, these can be immediately trans- 
lated into the winding numbers.) There are two immediate 
consequences of this conservation law. The first one concerns 
the ergodicity of the quantum dynamics: even for arbitrary 
local ring exchanges, the configurations corresponding to the 
different topological sectors (i.e., different overall slopes) can- 
not be mixed. Secondly, certain configurations, namely those 
corresponding to a maximal possible slope in at least one di- 
rection must be frozen since any local change in the height 
field would result in exceeding the maximal possible value of 
the step somewhere around the affected plaquette(s). 
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FIG. 2: Recipe of how to map dimer and FPL coverings to a height 
model. First the links of the bipartite square lattice are oriented in 
such a way that the arrows are always pointing from the black to the 
white sublattice. Next we start from one arbitrary plaquette and set its 
value to zero. We then set the height value of the other plaquettes in 
the following ways, (a) Dimer model: If the traversed link is empty, 
we change the height by Ah = 1 if the arrow is from left to right 
and by A/i = —1 if the arrow is from right to left. If the traversed 
link is occupied by a dimer, we change the height by Aft = —3 if the 
arrow is from left to right and by A/i = 3 if the arrow is from right to 
left. The dashed line shows one possible path assigning the heights to 
the plaquettes (b) FPL model: The rules are nearly identical to those 
for dimers except that | Afe| = 1 for both occupied and empty links. 
Since an FPL state is always an overlap of two dimer states, one can 
also obtain the same heights by simply adding those corresponding 
to the dimer states and dividing the result by two. 



For the fermionic variants of the FPL and dimer mod- 
els, there are two important differences which result from 
fermionic statistics. Firstly, only ring exchanges involving 
odd numbers of particles are allowed (as discussed in Sec- 
tion In particular, the smallest resonating plaquette has 
perimeter six, rather than four, as captured by Eqs. (O and 
and thus the height updates involve two neighboring plaque- 
ttes rather than a single one. From the point of view of the fun- 
damental conserved quantity of quantum height models - the 
global slope Hx{y) and Ky{x) of the height field - this differ- 
ence is immaterial since the overall slope cannot be changed 
by any local dynamics. In addition there is an accidental sym- 
metry in our model, specific to the quantum dynamics in our 
case: it separately conserves the total numbers of vertical and 
horizontal dimers (moreover, it conserves the number of hor- 
izontal dimers in even and odd rows separately, and similarly 
for the vertical dimers). Hence, all states can be classified with 
respect to these quantum numbers which are different from the 
aforementioned conserved slopes. Notice that these additional 
conservation laws will not hold if higher-order ring-exchange 
terms are added to the effective Hamiltonian. 

Secondly, the sign of the ring-hopping term in the fermionic 
FPL model depends on the occupancy of the middle bond. 
Both differences have important implications: the first one - 
for the ergodicity of the quantum dynamics, while the second 
one leads to a sign problem in the quantum Monte Carlo dy- 
namics. As we will see in the following subsection, the sign 
problem can be cured in certain cases. We note that the overall 
sign of iJeff (i.e., the sign of g) in Eqs. (O and Q is a mat- 
ter of convention (as has been first pointed out in Ref. S for 
the related bosonic QDM). Introducing a factor of i for each 
vertical dimer in an odd row of the square lattice and for each 



horizontal dimer in an odd column changes this sign. 
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B. Sign problem: Topological approach 



(b) 



I I 



I I I I 
I I 



u. 



■ 1 - - 1- ■ 
I I 

u 



1^1 ' 



I I I 



(d) 



o. 



!□ iXi CI 



- tKi =1 



= iXi =1 



1 




C IXI =1 



FIG. 3: (a, b) Representation of configuration as fully-packed di- 
rected loops, (c) Action of a "non sign-changing" process of the 
effective Hamiltonian (d) Three different actions of "sign changing" 
processes on the topology of loop configurations. 

We will now address the sign problem manifested in the op- 
posite signs of the A A and B B terms in Eqs. ^ and 
show that it can be avoided in certain (but not all) cases. Let us 
consider the surface with an even number of fermions remain- 
ing fixed at the boundary ("fixed" boundary conditions). Then 
fermion configurations are represented by closed loops as well 
as arcs terminating at the boundary. With no fermions at the 
boundary, there are closed loops only, as shown in Figure[3la). 
We add orientations to the loops as follows: (i) color the areas 
separated by the loops white and grey, with white being the 
outmost color; (ii) orient all loops so that the white regions 
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are always to the right. In the presence of arcs. Figure [3lb), 
we can choose how to close them on the outside without inter- 
sections (the outside region has no dynamics). Letting white 
be the color at infinity, we orient the loops as described above. 

We now consider how the effective Hamiltonian affects the 
loops covering. We first notice that the B ^ B processes do 
not change the loop topology (or their number) while A ^ A 
flips always do, with the possibilities presented schematically 
in Figure[2c,d). The relative signs resulting from the A ^ A 
flips are consistent with multiplying each loop configuration 
by i\—iY, where r (I) is the number of the (counter-) clock- 
wise winding loop. Hence, by simultaneously changing the 
sign of the A ^ A flips and transforming the loop states 
\C) — i'''^^^ |£), we cure the sign problem thus mak- 

ing the system effectively bosonic. 
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olated can be seen from a counter example shown in Fig- 
ure Ub). Here the five fermions hopp around an odd number 
of fermions in the middle and thus yield a sign change. How- 
ever, the fifth order process here does not change the number 
or the orientation of the loop. Consequently, the sign gauging 
based on the topology of the loops breaks down once these 
terms are included. 



C. Sign problem: Algebraic approach 

So far, we have attacked the sign problem using topolog- 
ical and computational arguments. In order to complement 
these approaches with an algebraic expression of the gauge 
transformation, we are seeking a possible way of representing 
the Hamiltonian (O with the kinetic part in such way that the 
"sign problem", which is manifested in the relative sign of the 
two terms, disappears. 

We use grey-and-white coloring of the plane as has been 
discussed above. The two distinct possibilities of coloring a 
flippable rectangle (up to an overall color reversal) are shown 
in Figure|5] Notice that in both cases, both plaquettes forming 
a flippable rectangle (marked by the crosses) change their col- 
ors as a result of a "flip". The key observation though is that 
in one case the color of these plaquettes is different and in the 
other it is the same. 




FIG. 4: (a) GS energy and energies of lowest excited states of the 
effective Hamiltonian in sectors with different slopes (kx, 0). Left 
panel: Data from an exact diagonalization of J/cff on a 72 site cluster 
where the Fermi sign is taken into account. Right panel: Same data 
from a calculation with excluded Fermi signs, (b) A counter example 
of a fifth order process which changes the sign of the wave function 
but not the topology of the loops covering. 

Note that this construction need not work for the periodic 
boundary conditions. Firstly, only even- winding sectors on a 
torus allow two-coloring. Secondly, even for even windings, 
we can start from a given configuration and find sequences of 
hexagons-flipping processes which invert the coloring when 
returning to the initial configuration (not shown here). This 
means that it is possible to dynamically reverse the coloring 
while and thus the sign gauging fails for these cases. It ap- 
pears though that for the periodic boundary conditions on even 
tori (preserving the bipartiteness of the lattice), the lowest- 
energy states belong to the sector where such transformation 
works, as confirmed by the exact diagonalization results (Fig- 
urellta)). 

The presented non-local loop-orienting construction is fur- 
thermore restricted to the ring-exchange processes of length 
six. The next higher non-vanishing terms result from pro- 
cesses involving five fermions hopping around a polygon of 
length 10. The fact that the above derived rules for the sign 
change due to the processes involving three fermions are vi- 





FIG. 5: Two distinct "plaquette flips" and corresponding colorings 
of the plane. The two moves are different by the relative "— " sign. 
The dashed diagonal lines represent either of the two equivalent (for 
the purpose of our argument) locations of a dimer along one of the 
adjacent sides of a square. The red crosses indicate the the double 
plaquette around which the loop segments resonate. 

Let us define an SU(2) "isospin" cr defined on each plaque- 
tte so that the eigenvalue a' = ±1 indicates whether the pla- 
quette is white or grey. We can then rewrite the Hamiltonian 
^ in terms of these new variables: 



Hi 



kin 



(4) 



We will now discuss the terms and show that (|4|i is a good 
representation of the original effective Hamiltonian. The op- 
erators Aflip(fT^) is a projector which annihilates a state un- 
less a legal flippable configuration exists in the plaquette pair 
(i, j), in which case its eigenvalue is 1. The explicit form of 
Aflip(cr^) is complicated - but what is important, it can only 
depend on the ct^'s of the two affected plaquettes and their 
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neighbors. (This is just a formal consequence of the fact that 
if we can look at the coloring and tell a flippable configuration 
from a non-flippable, hence there should be an operator that 
can do the same thing by "measuring" such coloring which in 
turn is uniquely determined by the values of cr^). If Afiip(tT^) 
returns 1 upon acting on a configuration, then the operators 
afaj simply flips the two isospins (i.e. flips the colors of 
both plaquettes), i.e., it inverts the colors of the two plaque- 
ttes. The operator erf cr| determines the overall sign of this flip 
by distinguishing the two cases shown in Figure |5] In other 
words, it yields the sign change which occurs in the B pro- 
cesses and represent the sign problem that we are attempting 
to tackle! 

The resolution of the sign problem in the chosen represen- 
tation is actually trivial: a^a^ = ia^ and therefore 



irk.„ = -.9^fTf'T|Afl,p(a--). 



(5) 



But the only "physical" direction in the isospin space is z 
- it corresponds to the plaquette being grey or white. The 
choice of x and y directions is a "gauge" choice. Hence we 
can always do a rotation around the z axis, so that y ^ x {in 
fact, we should rotate y — > ±x for the sites of the even/odd 
sublattice in order to consume the extra minus sign in Eq. (|5])). 
The result is 



(6) 



which is precisely Eq. ^ but without the different signs! No- 
tice that this "spin rotation" affects the phases of the wave 
functions (written in ct^ representation). This construction 
coincides with the "loop transformation" we have presented 
earlier. 

The non-local nature of this transformation is obvious from 
the fact that cr's are non-local operators: 



exp 



•E 

fcec 



(7) 



where C is a cut going from the given plaquette to the bound- 
ary (assumed white for simplicity), and Ck are creation 
and annihilation operators for fermions at wavenumber k. 
Sfeec '^\^k simply counts the number of dimers crossing this 
cut (i.e., the number of fermions on the cut). This is actu- 
ally identical to the Jordan-Wigner transformation in ID! This 
suggests that our model is "secretly" ID-like, which explains 
why the fermionic sign can be gauged away. Notice that this 
nice property appears to be a consequence of restrictive three- 
particle dynamics described by Eq. (O; adding five-particle 
ring exchanges seems to lead to a genuine sign problem that 
cannot be similarly fixed by a gauge transformation. 



IV. ROKHSAR-KIVELSON (RK) POINT 

Numerous examples demonstrate that a Hamiltonian con- 
taining only kinetic, plaquette flip (or ring exchange) terms 



leads to an ordered phase stabilized by the "order by disor- 
der" mechanism. On a square lattice, both the usual QDM 
and the quantum 6-vertex (FPL) models are at their plaquette 
phas o~°-~''"*"'"*' characterized by a checkerboard pattern (half 
as dense for the QDM case) of resonating plaquettes which 
break translational but not rotational symmetry of the lat- 
tice. A vary similar plaquette phase occurs on a honey- 
comb lattice—, where the dimer model and FPL models are 
simply dual to one another Moreover, a somewhat more 
complicated broken symmetry state of the quantum dimers, 
dubbed -\/l2 x ^/V2, has been found on a non-bipartite, tri- 
angular latticeJ^i^ In the fermionic case at quarter and half- 
filling the Hamiltonian ^ has also been shown to have an 
ordered ground state and fractionally charged, "weakly" con- 
fined excitations. 

In order to explore the possible quantum phases with 
fermionic ring exchanges, we follow the example of Ref. 
and extend our model's parameter space by adding a potential 
term as follows: 



h: 



'(1/4) 



eff 



E 



-g{\A){A\ + \A){A\) 



^^ji{\A){A\ + \A){A\)] (8) 



for the quarter-filled case and 



^eff = E [~9 {\A){A\ + \A){A\ + \B){B\ + \B){B\ 

V 

+/i {\A){A\ + \A)^\ + \B){B\ + (9) 

for the half filled case. Here we use the short-hand notations 
introduced in Eqs. Q and (O and in both cases applied the 
gauge transformations discussed in Section lTlIBI to change the 
sign of coupling constant in front of the A -i^ A terms. The 
Rokhsar-Kivelson (RK) point corresponds to fi — g where the 
Hamiltonian can be written as a sum of projectors 



(10) 



H^^\,=,=gY,[[\A)-\A))[{A\-{A\ 
p 

+ {\B)-\B))[{B\-(B\ 



for the half-filled case; only the the first line, i.e., A terms 
are present for the quarter-filling. The important feature of 
this Hamiltonian is that it is positive semidefinite, hence (i) 
a zero-energy state is automatically a ground state and (ii) 
such a state must be annihilated simultaneously by all projec- 
tors. The consequence of the fact that the gauge-transformed 
Hamiltonian is of the Frobenius type is that for each sub- 
space of states connected by the quantum dynamics, the equal 
amplitude superposition of such states is a unique "liquid" 
ground state. There is, however, a separate class of additional 
ground states - those not connected to any other states. These 
states by definition have no flippable configurations and hence 
are trivially annihilated by Hamiltonian ( flOl i. We note that 
such "frozen" states remain zero-energy states of Hamiltoni- 
ans (Is]),® for all values of fi, and remain the ground states for 
all II > g as these Hamiltonians retain their positive semidef- 
inite nature in this region. 
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FIG. 6: Fragments of possible "frozen" states: (a,b) are fluctuation- 
less under any local ring exchanges at quarter-filling (a) and half- 
filling (b); (c,d) are maximally flippable for bosonic, but fluctuation- 
less for fermionic ring exchanges; (e,f) are the examples of generic 
fluctuationless configurations of fermions at both quarter-filling (e) 
and half-filling (f); (g,h) are other examples of fluctuationless con- 
figurations of fermions at half-filling. 



While the fermionic nature of our original Hamiltonian has 
been seemingly gauged away, leaving us with a Frobenius 
type, we should recall that it is still manifest in the fact that 
only odd ring exchanges are allowed. This has a profound 
consequence for frozen states. All bosonic frozen states (i.e., 
frozen under a single plaquette flip dynamics) remain frozen 
in our case. This is immediately obvious on a superficial level: 
all rectangles flippable by either A A or B ^ B terms 
contain a flippable plaquette, hence a state with no flippable 
plaquettes is impervious the dynamics of Hamiltonians ^ or 



(|9|l. On a deeper level, bosonic frozen states remain frozen 
under fermionic dynamics since they are unaffected by any lo- 
cal ring exchange (see Figures|6|a,b)). In the height language, 
these states maximize the slope in some direction, making any 
local dynamics impossible. There are typically ~ 2^ such 
configurations, their exact number may depend on the bound- 
ary conditions. Such counting is easiest at half-filling where 
one may maximize the slope of an i x i system in in either 
zLx or ±y direction, which still leaves one with a freedom of 
2^ possible arrangements of steps in the other direction result- 
ing in a 4 X 2^ number of frozen states. In the quarter-filled 
(QDM) case there are four staggered states of the type shown 
in Figure |6(a)| that maximize the slope in either ±x or ±y 
direction, with the slope in the other direction automatically 
fixed to zero. Additional frozen states can be obtained by cre- 
ating diagonal domain walls between horizontal and vertical 
staggered states. Once again, all these states remain frozen 
under any local dynamics, fermionic or bosonic. 

The non-trivial consequence of the fermionic dynamics is 
the appearance of a large number of additional frozen states 
that is not related to maximizing the slope, in fact many of 
these states are in the flat sector. Two such examples are 
presented in see Figures |6lc,d). Both the columnar state of 
quantum dimers (Figure [6(c)]) and the analogous state of FPL 
(Figure [6(d)| i actually maximize the number of flippable pla- 
quettes, yet both are fluctuationless under any fermionic dy- 
namics. These states (along with staggered states) are mem- 
bers of a larger family of frozen fermionic states shown in 
Figures |6te,f). At quarter-filling, any state with only horizon- 
tal or only vertical dimers clearly cannot be changed by the 
application of the | A)(j4| terms since such terms necessarily 
involve both horizontal and vertical dimers. Moreover, it is 
easy to show that any ring exchange consisting of alternat- 
ing occupied and empty bonds will involve an even number 
of fermions in this state and hence is not allowed. There are 
2x2^ such states. A similar logic dictates that all FPL states 
like that shown in Figure [6(0| are frozen under any fermionic 
dynamics: whenever the loops can be sorted into two sets 
(shown in red and blue) in such a way that an empty bond 
always connects loops belonging to different sets, all ring ex- 
changes are necessarily even. The number of such states has 
been counted by us in Ref. |0] and shown to scale as 4^"*"^ for 
large L. A number of additional fluctuationless states with the 
same "bipartite" property of loops - see e.g. Figures |6jg,h) 
- can be constructed, yet we were unable to systematically 
count them. Hence it remains an open question whether our 
RK point is characterized by a t ruly extensive GS entropy like 
that found in Ref. IH or Refs. lUiS. 

In particular, one wonders whether our model can be actu- 
ally related to this latter class of supersymmetric models2ii2i 
in any simple way. In short, the supersymmetric Hamiltonians 
are constructed as 



HsvsY^{Q + Q^f = {Q,Q^} 



(11) 



where — ~ 0. It has been shown2ii2i that 

if one chooses Q — where P^j^ projects out 

the states with fermions occupying sites adjacent to i, then 
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such Hamiltonians may have an exponentially large num- 
ber of zero-energy states (which are necessarily ground 
states). This conclusion, however, is sensitive to both 
the filling fractions and the underlying lattices. While 
the checkerboard lattice has actually been studied in this 
context^^ - under the name of "square dimer lattice" - our 
Hamiltonian is clearly different. The ring exchange terms 
can in fact be produced within this framework by choos- 
ing Q = E«ij,fe,/,m,„» {CiCkCmPjPlPn ± CjClCnP^PkPm) 

where Pi ~ 1 — Hi; the sum is performed over all "flippable" 
rings consisting of sites i,j,k,l,m,n. The resulting super- 
symmetric Hamiltonian is both local and contains ring ex- 
changes, yet it also contains other terms (including additional 
one- and two-particle hopping terms). What is the relation be- 
tween the states of such a Hamiltonian and our Hamiltonian at 
its RK point is an interesting open question. A related ques- 
tion is whether our Hamiltonian can actually be written in a 
supersymmetric form by a more judicious choice of Q. 

Whether or not this is the case, it remains quite suggestive 
that the explanation for the large number of additional frozen 
ground states lies in a more restrictive nature of the fermionic 
dynamics. We reiterate that these "new" frozen states are not 
related to any tilt constraints in the height language. The fact 
the frozen states now include the maximally flat states (Fig- 
ure |6|c,d)) is particularly interesting since they are supposed 
to dominate the ordered phase of both quarter-filled and the 
half-filled bosonic model for ji < ji^hy maximizing the num- 
ber of the flippable plaquettes. At quarter-filling, this is a 
columnar phase quantum dimers; /ic ~ O.65 numerically/' 
At half-filling, this is a conventional anti-ferroelectric phase 
of the six-vertex model referred to as the Neel phase in 
Ref. IMl and as DDW phase in Refs. lfl9ll2l|] . its upper bound- 
ary has been found to be /ic ~ —0.3743122 The aforemen- 
tioned discussion makes it doubtful that such phases exist in 
the fermionic model. At half-filling, the phase dominated by 
a maximally-flat configuration shown in Figure |6(d)| appears 
to be replaced by the somewhat analogous "squiggle" phase 
shown in Figure [7(b)l '*^ In the quarter-filled case, the columnar 
phase appears to be replaced by a two dimensional analogue 
of the "R-state"i^^^ shown in Figure [7(a)] 




(a) (b) 



FIG. 7: Fragments of states stabilized by quantum fluctuations: 
(a) The state maximizing the number of three-fermion ring ex- 
changes at quarter-filling; (b) The "squiggle" state that maximizes 
the number of three-fermion ring exchanges at half-filling. The flip- 
pable rectangles involving the central plaquette are shaded. 



A. Excitations at the RK point 

The quantum dynamics of the Hamiltonian iJeff and simi- 
larly iJgff can be described in terms of a height model as de- 
scribed above.— The associated conserved quantities and 
Ky, as well as their insensitivity to a constant shift of the height 
field imply gapless hydrodynamic modes with a;(fc) ^ fc^ as 
fc 0. The liquid state of the FPL model at the RK point 
corresponds to the rough phase of the height model, in which 
case the modes can be identified as the so-called resonons at 
wave vector (tt, tt)^ and the equivalent of the piOns, here at 
Q = (0, 0).i^ilM^ In the following we will elaborate on this 
physics. 

For constructing gapless excitations in our model, we adopt 
the single-mode approximation^" which has been successfully 
used for the QDM on the square lattice.--™ A trial state with 
a momentum which differs from that of the ground state by 
q is constructed. The momentum is a good quantum num- 
ber in a translationally invariant system. Then the variational 
principle employed on the states at that momentum yields the 
basic result that the energy of the trial state provides an upper 
bound on the excitation energy at that momentum. There- 
fore, as a matter of principle, one can use the SMA to prove 
gaplessness. In order to demonstrate the presence of a gap, a 
different method is needed. 

Let |0) be one of the aforementioned "liquid" GS at the RK 

point. The operator d^^~\Y) creates (annihilates) a dimer at 
position r in direction (polarisation) f. The density opera- 
tor rif (r) = {Y)d^ {y) has the Fourier transform rif (q) ~ 
"-t"!^) cxp(iq-r). We use the operators rif (q) to construct 
the states |q, f) = nf(q)|0) which are for q 7^ orthogonal 
to |0) . By using the variational method, the excitation energies 
i?(q, f ) are bounded by 

£^(q,T)</(q)/^..(q), (12) 

where 

/(q) = (0|[n,(-q), [F^,, n,(q)]] |0) (13) 

is the oscillator strength and 

SfM) = (0|n^(-q)n^(q)|0) (14) 

is the structure factor, essentially the Fourier transformation 
of the dimer density correlation function. 

It is important that both the oscillator strength and the struc- 
ture factor can be evaluated as expectation values with respect 
to the ground state. The ground state correlations encode in- 
formation on the excitation spectrum and some of them can be 
captured by the variational method which leads to the above 
expression. SMA is useful then because there are situations 
in which gaplessness is present in a generic manner. The be- 
havior of /(q) near g — > can then be used to determine a 
bound on the dispersion of the soft excitations. A finite /(q) 
accompanied by a divergence of the structure factor can be 
used to infer gapless excitations. Such a soft mode is a classic 
signature of incipient order 

In order to calculate /(q), we observe that the density op- 
erators commute with the potential energy term, hence only 
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the kinetic energy term contributes. Let us first consider 
the case where f = x. By using the commutation relation 
[df, Uf] ~ Tdf repeatedly, we compute /(q) due to the res- 
onating terms. The result for the necessary commutators in 
order to compute the oscillator strength are (for a square lat- 
tice with lattice spacing a = 1) reads: 

[n,{-{q + k)), [-gT^/^,ns((Q + k))]] 

= AgT^^^ exp(iQ • R) sin^HQ + k) • y). (15) 

and 

[n£(-(Q + k)),[-gT^/^,n4(Q + k))]] 

= -4gT^/^ cxp(iQ • R)(l + cos((Q + k) • y)) 

X (-l + cos((Q + k) -x)) (16) 

where T^^^ and denote the resonant motion in the 

minimal ring exchange process with double plaquettes ori- 
ented in a column or in a row and with the middle dimer oc- 
cupied (A kind) or unoccupied (B kind). 

If we write k = q — Q and we consider q around three dif- 
ferent values of Q = (0, 0), (tt, tt), (0, tt) and therefore small 
k, then we can expand the results for the oscillator strengths. 
Then, to leading order we find that /(k) (k x f)^ (sim- 
ilar type of expression as the one for the RK modelS). Nev- 
ertheless, this result alone is not sufficient to give the whole 
information of the gapless excitations as the second necessary 
ingredient, the structure factor, can be also zero along the di- 
rection where k is zero. 

Therefore, let us now consider the structure factor S^x (q) 
both for the FPL and dimer models independently. We use the 
expression first given in Ref.|4^ Interestingly, the main result 
is that Sxx (q) 7^ except along the direction q = [qx, tt) 
where it vanishes with the exception of (tt, tt). At that point 
we have the gapless excitations termed as resonons by RK. At 
(0, tt) both /(q) and Sxx{^) are zero, but their ratio remains 
finite. In addition, for the FPL Sxx{^) shows no singulari- 
ties. In that respect, the FPL model differs from the hard- 
core dimer model (describing our model at quarter-filling) for 
which Sxx{^) diverges logarithmically at Q = (tt, 0) The 
difference is due to a slower algebraic decrease with distance 
of the correlation function for hardcore dimer covering. The 
correlation function for the dimer model at distance r reads: 

c^4^(-ir+^^^^ + (-i)^^ (17) 

while, for the FPL model it decays with a characteristic expo- 
nent which is higher than 3^. 

We have also verified the results for the structure factors of 
our model ^^^(q) by means of Monte Carlo simulations for 
both cases, see Figure[8ja) for the FPL model and Figure[8tb) 
for the case of hard-core dimers. The results remain identical 
(and symmetrical in the indices involved) for the orientation 
f = y. 



(a) 




FIG. 8: Structure factor Sa for the (a) FPL and (b) hard-core dimers 
on the square lattice. 



V. U(l) GAUGE THEORY 

The effective Hamiltonian (O conserves the number of par- 
ticles on each crisscrossed square, i.e., the number of loops 
or dimers touching each site is conserved. This conservation 
generates a local U{1) gauge invariance, as it is usually the 
case in models with local constraints ifisll . The gauge struc- 
ture suggests that we can gain further insight with respect to 
low-energy excitations by writing our model as a C/(l) lat- 
tice gauge theory. The usefulness of this approach has al- 
ready been shown for the quantum dimer model (QDM) on 
the square lattice' strongly correlated proton systems^S, as 
well as for three dimensional spin systems. The gauge the- 
ory can be derived for the dimer model as well as for the loop 
model in a very similar manner 



A. Dimer Model 

We start by defining an integer variable rij (x) for each link 
(x, X + ej) of the square lattice. The coordinates of a lattice 
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O CD 

X X.^^ X x.^^ 



FIG. 9: Labeling of links on flipable double-plaquettes at position x 
with unit vectors ei and 62. 



site are denoted by x, and 6^=1,2 are unit vectors along the 
axes as shown in Figure|9] The states |{?ij(x)}) span an en- 
larged Hilbert space which has integer numbers for the links 
instead only zero or one. We can consider the states \ {nj (x)}) 
as eigenstates of quantum rotor operators fij (x) with eigenval- 
ues nj(x). 

In order to express the effective Hamiltonian in terms of 
the nj{x.), we introduce phases 4>j{'x.) £ [0, tt) on the links 
which are canonical conjugate to nj{x.). Using the fact that 

±i0j (x) act as ladder operators, we can write 



exp 



^eff = f^E(^^W- 



1 



-2g cos 

{m,B} 



±(/. (.18) 



Here, the argument of the cosine term contains the sum over 
phases 4>j{x) with alternating signs around the polygons of 
perimeter six (double plaquettes). In the limit U /g — > 00, all 
"non-physical" states are projected out and Eq. ( fTSI ) is a faith- 
ful representation of the effective Hamiltonian Q. Next, we 
introduce staggered gauge and electric fields on the bipartite 
square lattice x = (xi, 2:2) by 

i,(x) = (-l)-^+-^0,(x) 

^,(x) = (^n,i^)^^y 

The local constraint that each site is touched by exactly one 
dimer reads 



= 



(19) 



(20) 



(a, 4 (x) - p(x)) iPhys. 
with the staggered background charge 

P(X) = + 

The lattice divergence is defined as 

AjEj (x) = ^i(x) - ^i(x - ei) + ^2(x) - ^2(x - 62). 

The dimer hardcore dimer constraint is now reflected by the 
standard Gauss' law ( [19] ) in the presence of a staggered back- 
ground charge density. Using the staggered fields as defined 
above, the Hamiltonian ( fTSb in reads 



^eff=C^E^'W-2.9E' 



.Q 



The argument of the cosine term denotes the oriented sum of 
staggered vector potentials Aj (x) around double plaquettes. 



B. Loop model 

Let us now consider the loop model which corresponds to 
the half filled checkerboard lattice. The matrix-elements of 
the effective Hamiltonian ^ have different signs depending 
on whether the link in the middle is occupied or not. Using 
the gauge transformation presented in section I III B I we can 
rewrite the Hamiltonian in a representation which has only 
negative matrix elements and redo the derivation presented 
above for the dimer model. The only difference between the 
dimer and loop model is the form of Gauss' law which reflects 
the loop constraint. Using the same notation as for the dimer 
model, the local constraint that each site is touched by exactly 
two dimer reads 



AjEj\Phy&.) = 0, 



(22) 



(21) 



Thus the corresponding Gauss' law has no background 
charges. 

Eq. (1211 has similarities with the Hamiltonian of the com- 
pact quantum electrodynamic (QED) in 2 + 1 dimensions in 
which the considered charges correspond to fractional charges 
of e/2j^ Polyakov showed for the compact QED in 2 + 1 di- 
mension that it has an unique and gapped ground state. Two 
charges are confined and the energy grows linearly with the 
distance between the two charges. Our model shows impor- 
tant differences. The fields Ej{x.) are half integers instead of 
integers and the constraint selects configurations with a back- 
ground charge p(x). This leads to a frustration which is re- 
flected by the macroscopic degeneracy of the classical ground 
states. Furthermore, the definition the flux in the cosine term 
differs. 

The formulation (l2Tl i provides an excellent starting point 
for further systematic investigations. The plaquette duality 
transformation allows to map the Hamiltonian to a height 
model and to use path integrals for a detailed study of the 
ground-state as well as low-energy excitations^^. 



VI. CONCLUSION 

In this study, we analyzed in detail a model of strongly cor- 
related spinless fermions on a checkerboard lattice both at 
half and quarter-filling by mapping it onto a quantum dimer 
and a quantum fully-packed loop models on the square lat- 
tice respectively. The lowest-order effective Hamiltonian is 
given by ring hopping processed of three fermions around 
hexagons. The symmetries and conservations of the model 
were studied in detail. We found a large number of frozen 
states which are not related to any height constraints and spe- 
cific for the fermionic case. Using a topological and an al- 
gebraic approach, we showed that the sign problem, induced 
by the fermion statistics in the half filled case, can be cured. 
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We explained and studied different aspects of this observation 
and showed that it does not remain true when further ring ex- 
change processes are also taken into account. By adding an 
extra term to the Hamiltonian, we fine tuned the system to the 
Rokhsar-Kivelson point where the ground states are known 
exactly. Using the proposed sign gauging and the single mode 
approximation, we argue that the Rokhsar-Kivelson point is 
a quantum critical point and the system crystalizes once it is 
slightly detuned. Moreover we mapped the system at both fill- 
ings to the corresponding gauge theories which dictates con- 
finement of fractionalized defects created by an addition (or 
subtraction) of fermions. 

Recent studies of other lattices such as kagome at 1/3 
fermionic filling concluded that the defects (in the form of 
fractional charges) are confined.-^ This is also the case in the 
present study and it is a consequence of the confined nature of 
the gauge theory as explained above. The fact that the lattice 
under consideration is bipartite as well as the imposed local 
constraints, lead to the fractionalization of the defects with 
protected value of the fractional charge (one half in we dope 
the system with an extra fermion). 

Finally, we should caution the reader who may be left with 
an impression that aside from the fact that only odd ring ex- 
changes are present, the fermionic nature of the model consid- 



ered here can be essentially "gauged away". Therefore it must 
be reiterated that the methods we have proposed are of some- 
what limited applicability. Firstly, they only work for special 
filling fractions. Their non-generic nature is obvious from the 
fact that for the quarter-filled case, a simple local gauge trans- 
formation would suffice to make the quantum dynamics of the 
Frobenius type, while for the half-filled case a far more com- 
plicated non-local transformation was required. Therefore, it 
would be extremely interesting to study this model away from 
these special fillings and try to identify the features that are 
generic. 
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Such cancelation of even ing exchanges is not a universal property 
but depends on the lattice: e.g., distorting a square lattice (where 
the dimers live) into a rhombic one will change the relative am- 
plitudes of clockwise and counterclockwise exchanges and so that 
their contributions will not cancel any longer. 
The mapping between FPL and six- vertex configurations is done 
by subdividing the square lattice into even and odd sublattices and 
orienting occupied bonds from odd to even sites, with the opposite 
orientation for vacant bonds. 



